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Abstract 

We compute the radial distributions of the vector and axial charge density as 
well as of the matter density in the heavy-light mesons on the lattice. The results 
for the lowest lying static heavy-light mesons and their first excited counterparts 
are obtained with Nf = 2 dynamical quarks of Wilson (clover) type and with the 
various improved static heavy-quark actions. From these distributions we were able 
to compute the corresponding charges and slopes of the form factors (i.e. (?" 2 ))- 
Those results are presented, and briefly discussed too. 
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1 Introduction 



Numerical simulations of QCD on the lattice are bound to provide us with information on 
nonperturbative dynamics among quarks and gluons. Not only do they allow us to compute 
the hadronic quantities with controllable systematic errors, they also give us access to 
studying the physics of confinement in more detail. In that respect the static limit of QCD 
is particularly attractive. Many studies of interaction between the two static (infinitely 
heavy) quarks lead to the results that helped us shape our understanding concerning the 
interquark potential (for a review see ref. [1] and references therein). The internal dynamics 
of heavy-light mesons in the static limit is far more complicated than that of the heavy- 
heavy ones, but it is even more interesting with regard to numerous phenomenological 
applications. For example, one can calculate spatial distributions of the electro-magnetic 
charge and of the light-quark densities, which have a very direct connection with the 
wave functions of quark models. Now, let us recall that the quark models are still a 
necessary tool in understanding hadronic spectra and various other properties of hadrons, 
including hybrid states. In fact, the lattice data show a striking agreement with the wave 
functions predicted by certain firmly established quark models, as we shall discuss in our 
next paper [2] (namely, certain models with Dirac equation and others formulated in the 
so-called Bakamjian-Thomas framework). On the other hand, this spatial distributions 
constitute a detailed testing ground for the quark models, for their choice of parameters, 
as well as for the basic assumptions, and therefore may be a precious means to select among 
them and/or to improve them. 

Recent progress in lattice QCD makes it possible to tackle this physics in more detail. 
In this paper we will compute the radial distribution of the matter, of the electric charge 
and of the axial charge by calculating the suitable hadronic matrix elements, which allow us 
to study these distributions as the light-quark is pulled away from the static source of color. 
Of essential importance is to have a good behavior of the static heavy-quark propagator, 
i.e. that the statistical noise does not become overwhelming as the separation between 
the source operators is increased to moderately large times. In recent years, substantial 
progress in that respect has been made. We will discuss the implementation of various 
improved static heavy-quark actions on the lattice, and comment on the benefit of each 
one of them. In producing the heavy-light mesons we use the light-quarks with mass that 
is varied between 0.7 x mf ys and 2 x mP hys , with mf ys being the physical strange quark 
mass. The gauge field configurations produced by the CP-PACS Collaboration [3] include 
the effect of Nf = 2 dynamical quarks, the masses of which are equal to that of the valence 
quark whose dynamics we will be studying in this paper. Besides the usual ground states, 
i.e. heavy meson doublet with the quantum numbers j p = (l/2)~, we will also study the 
similar properties of their orbitally excited counterparts, j p = (l/2) + . 

In the course of the following presentation we will occasionally make reference to some 
phenomenologically interesting information that comes out from our computation. How- 
ever, a more detailed discussion concerning the comparison of our results with various 
quark models will be addressed in a separate paper. 

In Sec. [2j we give a brief description of the lattices that we use in this work. We then 
discuss various heavy (static) quark actions and discuss the results from the two-point 
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heavy-light correlation functions. Description of the calculation of the radial distribution, 
and the presentation of our results is made in Sec. [31 There, we also present our results 
for the integrated distributions, and we finally conclude in Sec. HI 

2 Heavy-light mesons: two-point functions 

In this section we will discuss various improvement schemes of the static heavy-quark 
action. We will then test the gain of such an improvement on the two-point correlation 
functions involving the heavy-light meson source operators. Before we embark onto that 
discussion we will first give a short description of the lattices that are used in this work. 

2.1 About the lattices used in this work 

We use the gauge field configurations produced by the CP-PACS Collaboration, with Nf = 
2 mass-degenerate dynamical quarks: more specifically, they are obtained by using the 
Iwasaki gauge field action and the 0(a) improved Wilson quark action [3]. From the 
available data-samples we use the ones that are obtained on the 24 3 x 48 lattices at (3 = 2.1 
for which the lattice spacing is a ~ 0.1 fm (or a -1 ~ 2 GeV). From the large sample 
of configurations we choose the ones that are fully decorrelated. We actually take the 
configurations that are separated by 80 HMC trajectories. As a first step, we reproduced 
the results obtained in ref. [3], for the fully unquenched combinations (i.e. those for which 
the mass of the sea and valence light-quarks are equal). The pseudoscalar masses and the 
corresponding decay constants are obtained from the fit to 

((?7o75?)o,o(?7o75?kt) — ► \ fp m P exp(-mpt) , (1) 

X 

at large time separations which are found for t e [10, 23], like in ref. [3]. In Tabled] we list 



^sea ^val 


r mp [ m P (GeV)] 


mp/my 


r fp/Z A [ f P /Z A (GeV)] 


0.1357 


2.425(11) [1.025(5)] 


0.802(4) 


0.649(8) [0.274(3)] 


0.1367 


2.101(9) [0.888(4)] 


0.751(5) 


0.599(8) [0.253(4)] 


0.1374 


1.800(10) [0.761(4)] 


0.690(8) 


0.576(17) [0.243(7)] 


0.1382 


1.332(10) [0.563(4)] 


0.472(12) 


0.525(18) [0.222(8)] 



Table 1: The values for the light pseudoscalar masses and decay constants obtained from the standard 
exponential fit to the two-point correlation functions, on the configurations produced by CP-PACS, at 
P = 2.1 with Iwasaki gauge and the Wilson clover (csw = 1-47) action. 

those results together with the corresponding vector meson masses, which we computed 
from the correlation function similar to that in eq. ([IT), after replacing 7075 — > 7$. We used 
the values of r^ja from ref. [3] and tq = 0.467 fm, to convert to physical units. 
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2.2 Static heavy-quark on the lattice 

The notorious problem that troubled lattice practitioners for many years was the fact 
that the signal of the heavy-light correlation functions with infinitely heavy-quark was 
prohibitively noisy even at moderately large time separations. To reach the ground state for 
a static heavy-light meson it was necessary to implement some kind of smearing and hope 
that the lowest lying state was indeed the one observed from the effective mass plots at short 
time separations between the smeared sources in two-point correlation functions. Checking 
whether or not the smearing procedure is efficacious, by confronting to the correlation 
functions computed with the local sources, was practically impossible. 

Over the past several years that problem has been overcome by various recipes to 
improve the static quark propagator, i.e. the Wilson line. The implementation of the 
so-called "FAT-link" procedure and even more optimized versions with hyper-cubic 
blocking (HYP) [5| resulted in an exponential improvement of the signal to noise ratio of 
the static-light correlation functions [6] , making the corresponding physics results far more 
reliable than they were in the past. In the following we briefly describe these procedures 
and illustrate the benefit in using them. 

The latticized version of the static heavy-quark effective theory Lagrangian (HQET), 
in the frame in which the heavy-quark is at rest, reads [7j 

£hqet = J2{ hf W [h{x)-V (x-b)ih(x-6)] + Sm(g 2 )h^(x)h(x)Y (2) 

where V^ix) is the lattice link variable, and we include the power divergent term of the 
discretized HQET [5m(g 2 )h^(x)h(x)} which is a regularization artifact. For notational 
simplicity we omit it hereto, but we will of course discuss its impact when appropriate. 
The static heavy-quark propagator is then 

s h (xAy,t) = (n *tf( n )) 5 (*-® ( 3 ) 

\n=0 / 

with Vo(n) = Vq(x, n). Concerning the Wilson line on the lattice we will distinguish the 
following cases: 

• The Eichten-Hill (EH) static quark propagator is the simplest form obtained by re- 
placing [7] 

V (x) - U* H (x) = U (x) , (4) 
with Uq(x) being the simplest (unimproved) temporal link on the lattice. 

• The FAT-link static quark propagator is the first improvement which involves neigh- 
boring staples and is obtained by averaging over the six staples that surround the 
link g] 

V (x) - f/ FAT - 6 (x) = \ £ Uf^(x, x + i). (5) 
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By that modification the Wilson line becomes a thin flux tube. All the properties 
of the corresponding Wilson line remain the same so that the systematics of lattice 
calculations in principle is not changed. 

The HYP static quark propagator correspond to a more efficient choice of the lattice 
static quark action, in terms of improving the quality of signal for the heavy-quark 
propagator, and is obtained by an optimized combination of links belonging to the hy- 
percube surrounding a given link. This so-called HYP procedure can be summarized 
by the following three steps: 



V (x) - t/ HYP (x) = Proj su(3) 
where the links V^ v {x) are defined via 



1 - a x )U Q {x) + yY1 V MVo;i(x + i)vU x + 0) 



V^ v (x) = Proj SU (3) 
and finally V fl . U)P (x) is defined via 



(1 - a 2 )U li (x) + ^ V P*A X ) V ~w( X + P) V p;vA X + V) 



V KV , p {x) = Proj SU (3) 



l-a 3 )Up(x) + ^ Yl U a (x)U„(x + a)Ul(x + n) 



±a^p,vp 



• (6) 



The optimal choice found in ref. [5|, in which the HYP procedure was actually pro- 
posed, is a = (a\, 02,013) = (0.75,0.6,0.3), and we will refer to it as HYP-1. The 
other optimal choice, a = (1.0, 1.0, 0.5), was proposed in ref. [8] after minimizing the 
noise to signal ratio. We will refer to that option as HYP-2. 

The choices for parameters a made in refs. [3 [8] were obtained from the numerical stud- 
ies with the Wilson plaquette lattice-gauge action. In this paper we use the gauge field 
configurations generated with the Iwasaki gauge action. Instead of looking for yet another 
optimal choice for the coefficients a, we will explore the above four possibilities, study the 
impact of each of the improved action on the static-light quantities, and then choose the 
one that gives us the best results in terms of the quality of the signal and the associated 
statistical errors. Only then we will be ready to study the three-point functions. To illus- 
trate the benefit of the above-mentioned improvement we compute the propagator of the 
lowest lying static-light mesons, i.e. the correlation functions 

C 2 {t) = Y (h{x)^q(x) {h{y)^q{y)) ] ^ = (h( x H<l( x ) {KvhiQ(y)) 

y i,y 

= 7n ( Tr hsSh(x,0] y,t)j 5 Sg(y,t;x,0)}) + (Tr [75^ (z, 0; y, -t)j 5 S q (y, -t; x, 0)] \,{7) 



where one source operator has been fixed at some x at t = and the other is free and runs 
over all the lattice indices. The first line in the above equation exhibits the heavy-quark 
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H 




EH 


FAT-6 


HYP-1 


HYP-1 2 


HYP-2 


HYP-2 2 


0.1357 


?S 2 


104(27) 


106(9) 


111(7) 


111(5) 


111(5) 


103(4) 




£9 


0.770(36) 


0.613(11) 


0.578(8) 


0.534(6) 


0.548(6) 


0.506(5) 




2 

^9 


49(29) 


59(14) 


63(13) 


62(9) 


65(10) 


55(7) 






1.07(10) 


0.911(40) 


0.877(34) 


0.824(26) 


0.839(26) 


0.786(22) 






0.285(90) 


0.282(38) 


0.282(33) 


0.280(25) 


0.280(25) 


0.275(20) 






0.341(50) 


0.256(12) 


0.257(9) 


0.234(6) 


0.215(5) 


0.213(4) 


0.1367 


?S 2 


87(28) 


110(7) 


108(6) 


106(4) 


107(5) 


94(3) 




£* 


0.730(44) 


0.586(10) 


0.543(8) 


0.497(6) 


0.511(6) 


0.465(5) 




?s 2 
^9 


36(27) 


57(11) 


62(10) 


58(7) 


58(6) 


52(6) 






0.96(12) 


0.838(32) 


0.806(26) 


0.745(20) 


0.751(19) 


0.706(19) 




\ 


0.166(90) 


0.241(26) 


0.249(22) 


0.238(17) 


0.233(17) 


0.230(16) 




z L 


0.303(54) 


0.241(9) 


0.236(7) 


0.214(5) 


0.199(5) 


0.193(4) 


0.1374 


75 2 


85(28) 


115(10) 


112(7) 


101(5) 


100(5) 


88(3) 






0.714(47) 


0.582(12) 


0.533(8) 


0.475(6) 


0.486(7) 


0.436(5) 




?S 2 
^9 


36(44) 


77(18) 


71(14) 


65(8) 


63(8) 


59(6) 




^9 


1.28(20) 


0.833(38) 


0.772(30) 


0.709(20) 


0.713(21) 


0.676(16) 






0.477(181) 


0.234(33) 


0.227(27) 


0.223(18) 


0.214(18) 


0.229(14) 




Z L 
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0.294(55) 


0.246(12) 


0.236(8) 


0.205(5) 


0.188(5) 


0.181(4) 


0.1382 


9-5 2 


102(42) 


100(9) 


94(6) 


91(5) 


94(5) 


82(5) 






0.706(56) 


0.535(12) 


0.485(9) 


0.437(7) 


0.453(7) 


0.403(7) 




?S 2 
^9 


35(22) 


77(21) 


72(14) 


73(10) 


79(11) 


66(8) 




^9 


0.80(10) 


0.753(43) 


0.698(29) 


0.655(22) 


0.673(22) 


0.621(20) 




A, 


0.102(74) 


0.197(37) 


0.188(27) 


0.201(19) 


0.195(20) 


0.209(17) 




Z L 

9 


0.277(62) 


0.207(10) 


0.196(7) 


0.179(5) 


0.168(5) 


0.162(5) 



Table 2: List of our results obtained by using all six static heavy-quark actions : binding energies 
of the static-light j p = (l/2)~ mesons (£ q ), and of the j p = (l/2) + ones (£ q ), the orbital mass splitting 
(A q = £q — £ q ). We also list the values for the couplings to the smeared sources, Zq , Zq , as well as the bare 
matrix element of the local static-light axial current in heavy-quark effective theory, Zq = |(0|67o759|-Bg)|- 
All results are given in lattice units. 
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spin symmetry while in the second we average over the charge conjugated correlators. To 
illustrate the benefit of improving the static quark action we consider the effective mass 
(binding energy) defined as 



S cS = log 



C2{t) 

[C 2 (t + 1) 



and plotted in Fig. [U Even though the improved Wilson line allows the signal to exist at 
larger time separations, the statistical errors are still growing with time and a smearing 
procedure that would allow us to isolate the ground state at smaller time separations 
is desirable. It is particularly important for the studies of the three-point correlation 
functions, which will be discussed in the next section. We use the same smearing procedure 
as the one used in ref . [HI \TU\ , with the same set of parameters. In Table [2] we present the 
results of the fit to 

C 2 (t) = Z 1 ^ x exp (S g t) , 1,1' e (L, S), (9) 

where L and S stand for the local or smeared source operators. The fits are made in the 
interval 7 < t < 11 for all the light-quarks. The reader should not be surprised that the 
binding energies differ for different static heavy-quarks. This is because the residual mass 
term, Sm, in eq. ([2]) is different for different heavy-quark actions. In Fig. [1] we see that 
with the EH action the (smeared) correlation functions are prohibitively noisy as soon as 
one reaches t > 6. With improved actions the situation is much better and the signal is 
still present for t > 12. Besides lowest heavy-light mesons we also studied the correlation 
functions of operators coupling to the orbitally excited static-light mesons, i.e. those with 
j p = (l/2) + . They are obtained in a completely analogous way to what we just described, 
after replacing 75 — > I. Those correlators and the results obtained from them will be 
denoted by an extra tilde symbol. The fits to C 2 (t) are made in the interval 6 < t < 9. In 
addition, in that same interval we fit the ratio 

~ exp(A,t) , (10) 

where A q = S q — S q = Ax/2+ — A1/2-, which is free from the power divergent term com- 
ing from the HQET Lagrangian as written in eq. (j2J). In ref. [6] it was shown that the 
signal-noise ratio is exponentially improving with the improved Wilson lines. We found it 
interesting to check on the relative error of the signal as the time separation between the 
sources is increased. We therefore simply compute [5C 2 (t)}/C 2 (t), with 6C 2 (t) being the 
statistical error, and plot that quantity for various improved heavy-quark actions. This is 
shown in Fig. [TJ We indeed see the same phenomenon as in ref. [6] , namely, that the signal 
with improved Wilson lines remains stable longer than in the EH case, and its relative 
statistical error remains under 10% for quite large time separations. 

By "HYP-1 2 " and "HYP-2 2 " we denote the static quark obtained by reiterating the 
hyper-cubic blocking procedure (applying it twice). That procedure does not bring any 
gain if we consider the two-point static-light correlation functions alone. This can be 
viewed as an indication that the HYP procedure is well optimized. From the same figure 
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Figure 1 : In the upper plot we show the effective mass from the static-light two-point correlation function 
as obtained by using 6 different heavy-quark actions and the light-quark with K q = 0.1374, which is 
somewhat lighter but close to the physical strange quark mass. The lower plot shows the (exponential) 
growth of the statistical error of the two-point function as the time separation between the sources is 
increased. Notice that the error on the correlation function with HYP-2 heavy-quarks grows in the same 
way as with the HYP-1 2 one. 



and from our experience with three-point correlation functions, however, it seems that the 
HYP-2 procedure is somewhat better optimized and that for a better extraction of the 
signal with three-point correlation functions, the HYP-2 2 static quark action will be the 
most interesting for our purpose. As a side remark, we should mention that we tried to 
apply the HYP procedure thrice but we did not see any extra benefit in doing it. 

Before closing this section we should also comment on our procedure to compute the 
light-quark propagator. We do not use the "all-to-all" procedure [IT], which would certainly 
be very convenient, but instead we use the fact that we compute with the sources smeared 
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along each direction to at most -R max = 5 spatial points, so that in eq. (J7J) 

t-x 

l + 7o 



C 2 (t) = (Tr 



[VoinhzSqix^t; x,0)7 5 



n=0 



) 



we choose u x" to be anywhere within the smearing radius, i.e. we invert each of our light 
propagators 31 times [2x5 along each spatial axis plus with x at the lattice origin (0, 0, 0)]. 



3 Heavy-light mesons: three-point functions 

We now discuss the three-point functions. Since the heavy-quark is infinitely heavy, i.e. 
static, we can study the radial distribution of the light degrees of freedom within the 
bound state, which we choose to be either the ground state doublet J p = [0~, or its 
first orbital excitation J p = [0 + ,l + ]. In other words, we will compute the correlation 
functions by pulling the light-quark away from the static source of color by means of one 
of the local operators such as the scalar density (qlq), vector current (<77o<?) or the axial- 
vector current {q^i^^q)- In this way we will compute the matter, electric charge and the 
axial charge distributions. Such a study was first proposed in ref. [12], and in recent years 
in ref. [13]. 1 Preliminary results with an improved heavy-quark action (improved static 
Wilson line) have been reported in ref. [H]. In order to study the above-mentioned radial 
distribution we will study the following correlation functions: 



-fixed 



C T (-t sl ,t s2 ;r) = ^ (h{x)^q(x) q{z + r)T q(z + r) {h{y)^q{y)) ] 
= Yl ( Tr ^ Sh (y> x hsS q (x, z + r)Y S q (z + f, y)75]) z „ fixed (12) 



t s 2~l 



l+7o 



VqH j 5 S q (x,t s2 ;z + r,0) T S q (z+f,0;x,-t s i)j 5 



where one source of the heavy-light mesons is fixed to t y = —t s x, the operator whose 
distribution we want to compute is fixed at t z = 0, and t x = t s2 is free. In the practical 
calculations we also exchange the roles of t x and t y , and fix t x while letting t y free. For 
z = this would be the usual way we compute the correlation function. For z ^ 
and anywhere on the lattice the stochastic method of all-to-all propagators is needed. As 
we already mentioned we restrain our z to 31 points around the lattice origin. In other 
words, we make 31 inversions per light-quark q, per gauge field configuration. Finally, the 
components of the vector f can assume any of the values —L/2 < r XjV)Z < L/2. We should 
also add that in the above expression T = 70, 1. For large time separations in eq. (fl2l we 
then have 

Mrl = , S ^ 2 T{ Y;;/ I 7, > (B q \qTq(r)\Bq) , (13) 

(ZS) x exp[-Sg {t s2 - t sl )\ 



1 In fact the earlier study in ref. [12] was rather focused in extracting the gauge dependent Bethe-Salpeter 
amplitude, which -broadly speaking- also gives indication about the spatial distributions. 
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where we divided the couplings of the lowest states to the (smeared) source operators (2^) 
evaluated from the study of the two-point functions. One of the criteria to fix the correct 
fitting interval to extract the matrix element was to make sure that the result coincides with 
what is obtained if the three-point function is directly divided by the two-point functions. 
From that analysis and for t s2 — t s \ G {10, 11, 12, 13} we observe good plateaus of the ratio 
i?r(f), where one of the sources is fixed at \t s2 \ = 5 or at \t s2 \ =6 away from the origin 
t = 0, at which is placed our light-light-quark operator. For the sake of parity invariance, 
the above expression is exactly equal to zero if T = 7^75 (i = 1,2,3). Therefore, for the 
axial charge one of the sources should be the vector current, and we consider 



C yiTo (-t sl ,t s2 ;r) = ^ (h(x)~f 5 q(x) q(z + r)^a 5 q{z + r) (h(y)jiq(y)) 



! 

z— fixed 



x,y 



C (—t 1 t n-r) \t s i\,t s2 >o 

RwW = ,^2 Z m > (B q \qia^\BD , (14) 

[2$) x exp[-S q (t s2 -t 8l )\ 

where we obviously used the fact that the heavy-light vector and pseudoscalar mesons 
are degenerate. We considered all possible values of r 2 . We verified that the signal for 
various combinations of the spatial indices corresponding to the same r 2 are compatible 
among themselves, so that we could average over all of them. For larger values of r 2 the 
signals deteriorate and indeed after r = |f | ^ 12 the signal is poor and in some situations 
completely lost after r = |f | ~ 15, i.e. fully dominated by statistical errors. 2 This is not 
a worry because each of the distributions that we measure here are the rapidly decreasing 
functions, which eventually die out at larger values of r . Special attention should naturally 
be paid to lighter quarks for which the distributions are more spread out over the lattice. 
It appears that the quark masses that we deal with in this paper are not small enough to 
be sensitive to this kind of finite volume effects. 

The matrix elements defined in eqs. (TlBl) and (fl4j) are the scalar, vector, and axial 
densities, which we will denote by fr{r)- We will follow the terminology of ref. [13] and 
call / 7o (r) the charge density, /i(r) the matter density, and finally the one with f 7a6 (r) the 
axial charge density. 

As we already mentioned, here we also study the properties of the orbitally excited 
mesons, namely, the (l/2) + -doublet. To do so we repeat the same procedure as the one 
sketched above, by replacing in eq. (fT2l) 75 — > I, thus obtaining Cr(— t s i, t s2 ; r) {T = 70, 1), 
in addition to 

C 7i y 5 {-t sl ,t s2 ;r) = y2(h{x)lq(x) q{z + f)jij 5 q{z + f) (Hyhi^qiy))/ ,(15) 

— > ' '—fixed 



x,y 



where t z = and z = or any other of 30 points surrounding the lattice origin, as 
mentioned after eq. ()12j) . Desired distributions are then extracted from the ratios 

(2j?y x exp[-S q {t s2 - t sl )\ 



2 There are 7150 combinations of {rx,r y ,r z ) for \r\ < 12. After averaging over the equivalent combina- 
tions we end up with 122 points. 
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t s i\,t s2 >0 



> (B* 0q \q lil5 q(r)\B lg ) = f 7m (r) ,(16) 



{Zf) 2 x exp[-S q (t s2 - t s i)] 



where we used the standard spectroscopic notation according to which the states belonging 
to the (l/2) + -doublet, are labeled as Bq (scalar) and B\ (axial) states. We also keep track 
of the light-quark flavor when labeling the states. We should stress that the S-states in 
this work refer to the infinite 6-quark mass limit. As expected -from the discussion of 
the two-point correlation functions- the signal with the external excited states is less good 
than with the lowest states [(1/2) - ]. 

At this point we dispose of quite a number of distributions in r: we have two kinds of 
external sources [(l/2)~ and (l/2) + ] per light-quark mass, 4 light-quark masses, and for 
each of these we have 3 various densities that we are studying in this work. On the top of it 
we are computing everything with 5 different heavy-quark actions (FAT-6, HYP-1, HYP- 
l 2 , HYP-2 and HYP-2 2 ). That makes altogether 120 distributions with 122 points each, 
which would result in very cumbersome tables. 3 Before we opt for one of the heavy-quark 
actions that we consider in this work, we need to make several important observations: 

• In Fig. [2] we show the typical example of distributions f-y ,j ilSl i(r). Since they are 
rapidly varying functions it is more informative to multiply them by r 2 . Also, those 
distributions r 2 f^ 0tlil5 j(r) are shown in Fig. [2j Notice that the results are given in 
lattice units. We also explicitly show the value of each distribution at r = 0. The 
heavy-quark action used to produce the results shown in Fig. [2] is HYP-2 2 . Our 
observation is that by using various improved heavy-quark actions, the results are 
broadly compatible, but the quality of the signal with the HYP actions is much better 
than with the simple FAT-6 one. The resulting distributions obtained with the action 
HYP-1 2 are better than those obtained with HYP-1 (the shapes of distributions are 
smoother, and the statistical errors are somewhat smaller). The results obtained with 
HYP-2 action are as good as those obtained with HYP-1 2 (they are even numerically 
close). Furthermore, the results with HYP-2 2 are indistinguishable from the simple 
HYP-2 action. Therefore, it seems that with the Iwasaki gauge action the choice of 
parameters «hyp-2 = (1, 1, 0.5) is more optimized than the HYP-1 action. 

• The same observations are valid for the excited external states (l/2) + . Moreover, 
in this case the benefit of using the HYP-2 2 heavy-quark action is more obvious 
when applied to the data obtained with our two lightest quarks. For that reason in 
what follows we will restrain our attention to the results obtained by using the heavy 
HYP-2 2 action. 

• Before accounting for any renormalization factors we generically observe that for 
r = 0, /i(0) > / 707i75 (0). In addition, still for unrenormalized distributions, we 
confirm the observation made in ref. [13], namely that fi(r) is falling faster with r 
than / 7o (r). They cross at about r cross = 3 for our least light-quark, the point that 
then slowly shifts to r cross ~ 3.3, for our lightest quark. 

3 All those numbers can be obtained upon request from the authors. 
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• As we see from Fig. [2J all distributions r 2 /r(r) reach their maximum at some r m 
before asymptotically falling to zero at larger r's. From our data we see that, as the 
light-quark mass is lowered, the maximum, 

I [r7r(r)U= , (17) 

is reached between 2.4 and 2.9 for T = I, and 2.9 < r m < 3.7 for T = 7 . Finally for 
r = 7i75, we see that r m rs 3 for all our light-quark masses. 

For an easier use of our results, especially when comparing with various quark models, we fit 
the logarithm of our distributions to the polynomial of degree four, and the corresponding 
coefficients are given in Appendix A. 4 The explicit values for all distributions considered 
in this work, as obtained with the action HYP-2 2 , are given in Appendix B. 



3.1 Vector Charge 

In the following we explain how we compute the charge corresponding to each of the 
distributions discussed in this section so far. We checked that to a very good accuracy 
(indistinguishable by our error bars) 

r 2 f r (r)dr, (18) 

where on the left-hand side we have the sum over all of the lattice points. On the right- 
hand side, instead, to get the radial distributions we averaged over various combinations 
of {r x ,r y ,r z ) which correspond to a same value of r 2 . Above equality is not in principle 
obvious since the right-hand side makes use of the spherical symmetry of the problem, and 
we then integrate over the angular dependence, while our calculations are made on the 
cubic lattice. The reason for the validity of the above equality lies in the fact that our 
distributions are negligibly small even before they reach the edges of the lattice. Notice also 
that the upper bound of the above integral "oo" can be identified with r = 12 (in lattice 
units), and further we go beyond that point we are picking the contributions from the 
lattice points that do not respect the spherical symmetry. We tested that the integration 
by including the points that go beyond r = 12 never modify our integral by more than 1%. 
Finally, concerning the error estimate when doing the integral we estimate them by using 
the jackknife procedure. 

We now focus on the vector current and its charge. Since our heavy-quark is a spectator, 
the conservation of the light-quark vector current fixes the vector charges 

/•OO />OG 

# = 47t/ r 2 fl(r)dr, and = 4tt / r 2 /« o (r)dr, (19) 

Jo Jo 

to unity, once we multiply them by the appropriate renormalization constant, i.e. Zy (gfygv — 
Zv{go)gv = 1- This is therefore simply a method to nonperturbatively evaluate the vector 

4 The only exceptions are the scalar and axial distributions with the excited external states, which 
change the sign and therefore taking their logarithm would not make sense. 
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Figure 2: Radial distribution of the vector, axial and the scalar density within the static-light meson. 
On the left side the distributions are obtained with the external (l/2)~ mesons, and those on the right 
with the (l/2) + ones, light-quark is n q = 0.1374, which is around the physical strange quark mass. All 
results arc displayed in lattice units. 
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current renormalization constant Zv{gl), which we then confront to the perturbatively 
estimated values to chose the suitable "boosted" coupling. Only then we can use more 
reliably the perturbative formulas to renormalize other charges that will be discussed in 
what follows. In Table [3] we list our results for both g^p and g^ for all of our quark 
masses (i.e. K q ). Moreover, we present those values for all kinds of improvement of the 
heavy-quark propagator discussed in the previous section. We again see that the results 
with the action HYP-2 2 are indeed the best, and are compatible with those obtained by 
using other heavy-quark actions, as they should. It is also clear that the results with the 
excited heavy-light mesons (l/2) + are noisier, however compatible with what we obtain 
with the ground external states, (1/2) - . We also checked the saturation of the sum/integral 
in eq. (I19p . and we found that the integral is already fully saturated for r < 11, for all 
of the quark masses discussed in this paper. We illustrate this saturation for one of our 
light-quark masses in Fig. [31 To evaluate the renormalization constant we then linearly 
extrapolate l/g^ in am q to the chiral limit. 5 From our data with the HYP-2 2 action, 
the intercept is then Zy\g^) = 0.739(40). We checked, however, that this result remains 
stable regardless of the heavy-quark action in the spectator line. 

This last result in fact helps us choosing a convenient boosted coupling. It turns out 
that with the simple replacement g\ — > g 2 p = g 2 /(P), with the average plaquette being, 
(P) = 0.60220(6) [3J, we reproduce the nonperturbative result. Using the perturbative 
expressions reported in rcf. [16J, the coupling af (a) = 0.3775 and csw = 1-47, we obtain 
Z^\gl) = 0.780. From now on to renormalize each of our currents we will use the boosted 
perturbative values for the renormalization and improvement constants, the values of which 
will be reported when necessary. One interesting by-product of our analysis is the charge 
radius of the heavy-light meson. We remind the reader that the charge radius provides 
us with information concerning the slope of the electromagnetic form factor [F(g 2 )] of the 

5 Here am q = (1/2)[1/k 9 — l//c C rit]j where the critical hopping parameter has been evaluated in ref. [3] 
and reads K crit = 0.138984(13). 



11 



Heavy-quark 


K q 


<#> <r 2 >i- 




FAT- 6 


0.1357 
0.1367 
0.1374 
0.1382 


1.01(5) 17.1(2) 
1.08(5) 20.3(3) 
1.30(8) 23.2(3) 
1.24(7) 27.7(5) 


1.22(31) 17.3(18) 
1.12(22) 20.9(21) 
1.31(29) 17.2(25) 
1.43(36) 29.0(18) 


HYP-1 


0.1357 
0.1367 
0.1374 
0.1382 


1.05(4) 17.2(2) 
1.09(4) 20.5(3) 
1.26(5) 23.3(3) 
1.22(6) 28.2(4) 


1.23(26) 16.6(16) 
1.24(20) 21.0(17) 
1.24(23) 18.0(23) 
1.36(27) 29.4(17) 


HYP-1 2 


0.1357 
0.1367 
0.1374 
0.1382 


1.08(4) 17.7(2) 
1.13(3) 21.1(2) 
1.23(4) 23.8(3) 
1.26(5) 28.8(4) 


1.26(21) 17.7(13) 
1.20(16) 22.1(13) 
1.15(16) 18.0(20) 
1.43(23) 28.7(15) 


HYP-2 


0.1357 
0.1367 
0.1374 
0.1382 


1.08(4) 17.8(2) 
1.12(4) 21.2(2) 
1.21(4) 23.8(3) 
1.26(6) 28.8(4) 


1.24(19) 17.8(13) 
1.15(16) 22.7(13) 
1.11(14) 17.8(20) 
1.35(22) 28.7(16) 


HYP-2 2 


0.1357 
0.1367 
0.1374 
0.1382 


1.10(3) 18.4(2) 
1.15(3) 21.8(2) 
1.22(4) 24.4(3) 
1.29(6) 29.3(3) 


1.25(18) 18.4(12) 
1.19(14) 23.1(12) 
1.17(14) 18.9(17) 
1.44(22) 27.4(15) 



Table 3: Vector charges obtained by integrating the corresponding vector current densities and for both 
the (l/2)~ and (l/2) + external states. 



pseudoscalar meson P at q 2 = 0, i.e. 



(P(p>)\J™\P(p)) = F(q 2 ) (p + p') , <r 2 > = 6 



dF(q 2 



(20) 



q 2 =0 



where q = p — p', and J^ m = eQQ'y^Q + e^^^g, with eg^) being the electric charge of the 
heavy(light) quark. Similarly one can define the charge radius of the scalar meson. Since 
we already have radial distribution of the light-quark component of the electromagnetic 
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form factor we can simply compute its moment, namely, 



DC 



r fjo( r ) dr 

(r 2 )^ = ^ , (21) 

r2 Ao( r ) dr 





and the corresponding results obtained from our data are also listed in Table [31 Similarly 
for the (l/2) + -states. As it could have been anticipated, the saturation of the integral 
(lattice sum) in the numerator of eq. ( T2TT) occurs at higher values of r, since the extra r 2 in 
the sum enhances the importance of the values of the rapidly falling distribution / 7o (r) at 
larger radial distances. This is also illustrated in Fig. [3j We do not attempt an elaborate 
chiral extrapolation here. We relegate that discussion to our future publication. We note, 
however, that by simply linearly extrapolating to the chiral limit we obtain 

(r 2 ) t - = (31.6 ± 0.3) x a 2 =^ (r 2 )T- = 0.334(3) fm 2 , 

2 ? 

(r 2 ) i+ = (28.3 ± 1.6) x a 2 => (r 2 )T+ = 0.280(16) fm 2 , (22) 

2 2 

where we fixed the lattice spacing to a = 0.0995(4) fm, from r^ja = 4.695(18), as obtained 
on this same lattice [3], and by using rQ hys = 0.467 fm. It is interesting to note that these 
values are slightly lower but comparable to what one would obtain by applying the simple 
vector meson dominance, i.e. (r 2 ) VMD = 6/m 2 = 0.393 fm , and it is considerably smaller 
than the results obtained from the quark model of ref. [15]. 6 Another phenomenologically 
interesting result can be obtained by using our results ( 1221) to fix the left-hand side of the 
Cabibbo-Radicati sum rule for heavy hadrons discussed in ref. [T7] 



dF{q 2 



where the B* contribution has been transferred from the left-hand side to the sum on the 
right-hand side. By only keeping that very first term in the sum, we obtain a bound 

1 6 

r(£T± - B± 7 ) < -a em (r 2 ) \k,\ 3 = 1.16(1) keV . (24) 

ol 

Similarly, T(B*° —> B°j) < 0.29 keV, where we used |fc 7 | = (m%* — m 2 B )/(2m B *) = 
45.4 MeV. It should be emphasized that this upper bound is probably very weak as we will 
show in our forthcoming paper where we verify that, in the quark models, the magnetic 
moment is expected to be much smaller than what is suggested by this bound [2]. 



6 We would like to stress that this is a feature of that specific model. In our next paper we will show 
that other models are quite consistent with (r 2 ) i± obtained on the lattice [2]. 
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gf 


(r 2 )?- 

2 


gf 


(r 2 >?+ 

2 


gA 


(r 2 )t~ 

2 


g ( l ] (r 2 )t+ 

2 


0.1357 


0.938(33) 


12.6(3) 


0.91(17) 


26(3) 


0.695(18) 


15.1(2) 


-0.077(42) 


0.1367 


0.963(36) 


15.1(5) 


0.90(18) 


36(3) 


0.705(18) 


18.1(2) 


-0.035(40) 


0.1374 


0.954(43) 


15.3(7) 


0.70(22) 


34(6) 


0.683(18) 


20.7(3) 


-0.088(36) 


0.1382 


1.02(8) 


18.9(2.2) 


1.01(40) 


50(8) 


0.657(25) 


24.6(6) 


-0.220(55) 



Table 4: The scalar and the axial charges obtained by integrating the axial and matter densities as 
indicated in eqs. (|25I33|) . We also present the results for the corresponding mean square radius. The 
results for the external ground and the excited heavy-light states are displayed. The values presented in 
this Table refer to the static heavy spectator quark obtained by using the HYP-2 2 action. 



3.2 Scalar Charge 

Concerning the scalar charge and the mean square radius of the corresponding scalar form 
factor, we integrate the matter density 



4ir 



r 2 fj(r)dr, and ( r2 )f 



r 4 fi(r) dr 
r 2 fi{r) dr 



(25) 



and similarly for the (l/2) + states which we -as before- distinguish by an extra tilde. In 
this case we only present our results obtained with the improved heavy-quark action HYP- 
2 2 , and they are listed in Table HI The perturbative matching to the continuum includes 



the following constants: Z^ MS (1 / a) bpt 
extrapolating to the chiral limit we obtain 

0.606(37), and 

0.45(18), and 



0.602 and b h s pt 



1.472, 7 so that after linearly 



(r% 



= 0.193(9) fm 2 
0.54(6) fm 2 , 



(26) 



where we notice a strong deterioration of the signal for the (l/2) + external states. 

So far by "g" we denoted a light-quark without specifying its flavor. It stands for u or d, 
the two degenerate quarks in our configurations with iVf = 2 dynamical light-quarks. With 
the scalar density, however, one should be particularly careful. We combine our quarks 
in a flavor- doublet ip = [u d] T , so that the neutral current can be either flavor singlet or 
flavor triplet, i.e. 



ysing. 

J s 



xtrip. 
J S 



(27) 



7 As a side note we should state that for the situations considered in this paper, no 0(a) improvement of 
the bare current gives a non-zero contribution. The improvement of the renormalization constants includes 



Zs,v,a(9q) 



AO) 

J S,V,A 



(So 2 )[l 



l(0) 

^S,V,A\ 



(<?Q)am g ], with Z^y A {g$) being the renormalization constants in the 



chiral limit, while iif li(5o) are the linear mass counter-term coefficients. As we discuss in the text, these 
quantities are taken from boosted Tloop perturbative formulas [16) . 
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with T = diag(l, 1), and r 3 = diag(l, —1). In other words, 

Hu = Jfs- + jfP- , did = Jf g - - J| ip - . (2* 



it J 



The problem is that the computation with the singlet current contains disconnected di- 
agrams, which are particularly important in this case. Their computation on the lattice 
requires the implementation of the "all-to-all propagator" technique [IT] . In this work we 
neglect the disconnected diagrams (in other words, we compute only the flavor nonsinglet 
part of scalar density). Beside the computational complications in evaluating the discon- 
nected diagrams, an extra complication arise from the fact that, on the lattice, the qq 
operator mix with the cubically divergent (di/a 3 )H [and, away from the chiral limit, also 
with (c^wig/a 2 )!, and (^m^/a)!], where 0^1,2,3 are the constants whose values should be 
determined and mixing subtracted away before multiplicative renormalization. To check on 
the importance of disconnected diagrams we can compare our results with the light-quark 
mass dependence of the binding energy S q . In other words, we would like to verify if the 
following equality is satisfied: 8 

dS q _ /Bq] ^k ]Bq) = {Bq \qq\B q ) I g™(l/a) , (29) 



dm q dm L± 

where ri q = m u uu + m^dd, and either we choose q = u or q = d the evaluation of the 
left-hand side of the above equation will contain the disconnected diagrams. The last g$ 
on the right-hand side is our result obtained in eq. ff26l) in which the disconnected diagrams 
were omitted. If we fit our values for S q , obtained with the heavy-quark action HYP-2 2 
and presented in Table [2} as S q = So + pm q , we obtain 

S = 0.369(7) , p m (l/a) = 1.02(7) , (30) 

where we used the quark mass defined via the vector Ward identity, m q = [Z^ MS (1/ a)]~ l m q ° 
x(l — m q °^bs/2), and m q = 1/2(1/ K q — l//«crit)- This result is more than 60% larger than 
our value for g^ s (l/a) = 0.606(37). We also attempted using the quark mass defined via 
the axial Ward identity (the bare values of which are given in ref . [3] , which we recomputed 
and confirmed). The advantage of using the quark mass defined via the axial Ward identity 
lies in the fact that the correction due to bA~bp = —0.009, is negligible. From that fit we 
obtain that the slope is 1.28(9), and after accounting for the proper renormalization factor 
(i.e. [Zp(l/a)/Z A } hpt = 0.692), we have pi = 0.89(6), hence still about 50% larger than 
our g^ s (l/a) in eq. fT26l) . Therefore it appears that a contribution due to disconnected 
diagrams is indeed important and large. We checked that this conclusion remains as such 
if we fit around any of our directly accessed light-quark masses [i.e. S q = S\-\- p\{m q — mi)], 
since eq. f[2"9~j) is valid for any light-quark. 

The calculation of disconnected diagrams is not only demanding computationally. It 
also requires to solve a problem of mixing with the power divergent operator even in the 
chiral limit, which is therefore yet another source of systematic error when computing the 
light-quark mass dependence. As it is well known the binding energy computed on the 



3 This is also known as Feynman-Hcllmann theorem, discussed in ref. [18| and references therein. 
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lattice (S q ) cannot be directly compared with the continuum definition (A q ) unless the 
linearly divergent contribution due to the lattice regularization is properly subtracted, i.e. 
S q = A q + c((?o)/ a - in addition, when away from the chiral limit the subtraction constant 
becomes c(#q) = co((?o) + am q c i(9o)- By taking A q = A + p\ x Tn q , with A being the 
binding energy of the static-light meson in the chiral limit, then 

S q = A + ^-+[p A + c 1 (g 2 )]m q =► ^ = p A + c 1 (^ . (31) 

To check whether or not this latter artifact is (partly) responsible for not verifying the 
equality in eq. (1291) it is indispensable to compute the disconnected diagram and see to 
what extent its contribution would saturate the equality (129"!) . We should also add that 
eq. (I3"TT) is one reason more why -in HQET on the lattice- it is better to compute the 
mass differences, such as A q = £ q — S q because in that case the linear divergence cancels 
altogether. Linear fit to our data gives 

A q = 0.187(17) + 0.58(19) x mjp(2 GeV) . (32) 

In physical units this leads to Ao = 372(37) MeV. Notice however that the slope is (clearly) 
positive. Therefore, at the masses considered in this work , the lattice results exhibit a 
behavior opposite to what is established experimentally in the charm spectrum, namely 
that A q increases when one passes from the strange quark (q = s) closer to the chiral limit 
(q = u,d) (see discussions in [19] and references therein). 



3.3 Axial Charge 

In Table SJ we also presented our results for the axial charge and its corresponding mean 
square radius, defined analogously to what we did with the vector and scalar charges 

/*oo /»oo 

gf = 4tt / r 2 f A (r)dr , gf = 4tt / r 2 f A (r)dr , 

Jo Jo 

r 4 f A (r) dr / ^/^(r) dr 

, ( r % = -f^~: . ( 33 ) 

r 2 f A (r) dr / r 2 f A (r) dr 

Jo 

These quantities are particularly important because they are related to the P-wave soft 
pion emission/absorption in the transitions between the mesons belonging to the same 
doublet, as discussed in ref. [20]. In particular, 

(B q \q^ 5 q\B* q ) = gfe, (B* g \qj l5 q\B' lq ) = gfe* , (34) 

where e(e') is the polarization vector of the vector B* (axial B[ q ) static-light meson, and 
the axial current is now renormalized. After accounting for the axial renormalization 
Z A vt {g Q 1 ) = 0.803, and b A pt (gQ) = 1.346, the simple linear chiral extrapolation gives 

g A = g = 0.526(22) , g A = g = -0.157(48) . (35) 
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If this extrapolation is guided by expressions derived in chiral perturbation theory (HM- 
ChPT) |2I], we get 

SHMChPT = 0.444(18) , WchPT = -0.14(6) . (36) 

Furthermore, the simple linear extrapolation of the mean square radius to the chiral limit 
results in (r 2 )f_ = 0.246(3) fm 2 , and it measures the deviation of the axial coupling 

2 

to the pion that is pulled off its mass shell. In other words, we obtain that g(q 2 ) = 
g(0)[l + (r 2 )f_/6 x q 2 ], which is steeper than the shape predicted by the axial- vector 

2 

meson dominance, namely (r 2 )f_ = 6/m 2 ^ = 0.154(8) fm 2 . Since we are considering the 

2 

coupling to a pion (flavor nonsinglet) and not to the singlet states (77, 77'), the problem of 
disconnected diagrams is not present in this situation. 

We should also comment on the values of the axial coupling g^ , which straddle around 
zero and takes a more pronounced nonzero value as the light-quark mass becomes lighter. 
From the distributions that are presented in the Appendix A, and also the one already 
shown in Fig. El we can see that they change the sign and have zero at about r = 0.4 fm, 
so that in the sum the short distance (negative) and long distance (positive) contributions 
to the coupling cancel each other to a large extent, which is one of the reasons why 
the errors on g^ are much larger than in the (l/2)~-doublet case. 

4 Summary and conclusions 

In this paper we discussed various spatial (radial) distributions of the light degrees of 
freedom with respect to the static heavy-quark. We considered both the lowest lying heavy- 
light mesons, as well as the nearest orbitally excited ones. We used various improvement of 
the static heavy-quark action on the lattice. We observe the best efficiency achieved when 
the HYP-2 2 action is used. Concerning the (electromagnetic) charge distributions, we find 
that they have almost the same shape in both cases, i.e. when the external states are either 
jp = (l/2)~ or (l/2) + . The matter and the axial distributions have different shapes in two 
cases: while for the external (l/2)~ states they are strictly positive, for the (l/2) + states 
they change the sign: they start by being negative and change the sign at r pa 0.25 fm 
for the matter distribution, and at r pa 0.4 fm for the axial one. That zero moves slowly 
towards larger r, as we lower the light-quark mass. Such a qualitative difference of the 
internal structure of (l/2) + -mesons and the (l/2)~ states is particularly evident in the case 
of axial density. As a result, the integrated characteristics (axial charge) in (l/2)~ and in 
(l/2) + -states are completely different. These axial charges are also physically interesting 
because they describe the coupling of a soft pion to the (l/2)~ and to the (l/2) + states, 
respectively. They are essential parameters in the heavy meson chiral perturbation theory. 
From our analysis we found it to be ^n n . = 0.52(2), and ^HMChPT = 0.44(2), depending 
on the extrapolation formula used to go to the chiral limit. For the soft pion coupling to 
orbitally excited states we obtain, instead, gu n , = —0.18(5), and (ThmciiPT = —0.14(6). This 
further corroborates the fact that the model of chiral doublings of heavy-light hadrons [23] 
is so badly broken that even its qualitative features should be taken very cautiously. 
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For each of the distributions we computed the mean square radius, which is a new 
result. From the linear chiral extrapolation we obtain the values that are fully consistent 
with what one would naively obtain by using the nearest resonance approximation in the 
crossed channel. In particular by using our result for (r 2 }\- = 0.33 fm 2 , and the Cabibbo- 

2 

Radicati sum rule for the heavy-light mesons [T7] we obtain T(i?* ± — > B ±/ ~f) < 1.16(1) keV. 
Our results also indicate that the scalar charge and the matter density receive important 
contribution from the flavour singlet piece, namely, the one that contains disconnected 
diagrams, which have not been computed in this work. 

As a logical continuation of this project we plan to compare the distributions computed 
and discussed in this paper to those obtained by using two specific classes of quark mod- 
els [2]. An attempt to study the issue of disconnected diagrams, to complete our study 
of the matter density and make the comparison with dS q /dm q more conclusive, would 
be highly welcome. In this paper we did not extend our analysis to the radially excited 
states. Such an analysis may soon be possible thanks to recent improvement regarding the 
extraction of the properties of such states from the lattice [22J. A more elaborate study of 
chiral extrapolation requires the appropriate expressions in HMChPT, which we also plan 
to compute. Finally, even though our computations are made on relatively fine lattice and 
with O(a) improved Wilson quarks, it is necessary to check how if our result remain stable 
if even finer lattices are used. Such a study is underway. 

Phenomenologically more interesting is the comparison of results presented in this paper 
to various quark models. One such study is the topic of our forthcoming paper [2]. 
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Appendix A 



In this appendix we present the plots for all distributions fr{f)- Instead of plotting the results 
obtained with all 5 heavy-quark actions, we choose 3 of them: FAT-6 (Fig. Hj), HYP-1 (Fig. [5]) 
and HYP-2 2 (Fig. [6]). The reason is that the numerical results obtained with HYP-1 2 are 
very close to those obtained with HYP-2, and those with HYP-2 2 are better in quality than 
those with HYP-2 only in the situations in which the external state is orbitally excited, i.e. 
(l/2) + . Furthermore, to get a honest assessment of the quality of our results, we illustrate 
the distributions for our largest {n q = 0.1357) and our smallest {n q = 0.1382) light-quark 
mass. The plots are made for the distributions in lattice units and without accounting for the 
renormalization factors. To convert them to physical units, one should multiply the "x" axis 
by a = 0.0995(4) fm, and the "y" axis by 1/a 3 for fr(r), and by 1/a for r 2 f r (r). 

When fr{r) is strictly positive, we fit its logarithm to a polynomial of degree four, namely 
log[/r(r)] = Co + c\T + c 2 r 2 + C3?" 3 . Instead, in the cases in which fr{r) changes the sign, we fit 
to the forms listed below (the guiding principle is parsimony of parameters). 
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(co + cir 2 ) exp(-c 2 r 


1 15 1) 


(co + cir 3 - 5 ) exp (- C2 r^) 


0.1367 
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(c + c x r 2 ) exp(-c 2 r 


1.25 |) 


(c + Cir 3 - 5 ) exp(-c 2 r 


1.35 1) 



The coefficients from the fit within the interval r e [0, 12] for HYP-2 2 are listed below, and the 
fitted curves are also displayed in figs. 0], [5] and [6j 
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Figure 4: Distributions obtained with the heavy-quark actions FAT-6: upper set of two rows correspond 
to the external lowest states, and the lower set to othe excited states. In the upper row in each set the 
light-quark is our least light one (n q = 0.1357,), while in the lower it is our lightest (n q = 0.1382) . The 
black lines are obtained from the fits to the distributions fr(r). 
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Figure 6: The same as in Fig. [4| but with the heavy-quark actions HYP-2 2 
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Appendix B 



In this appendix we give the explicit values for the distributions /r(r) and fr(j~) which we 
call S, V and A for T = I, 70 and 7^75 respectively. We first present those obtained for the 
lowest heavy-light states, i.e. j p = (l/2)~, and then for the orbitally excited ones j p = (l/2) + . 
For each r we give simultaneously the values for all the light-quarks used in this work. 
We stress again that the light valence and the sea quarks are degenerate in mass. Finally, we 
need to explain how the tables are to be read. The first point (obtained at r = 0) is to be 
multiplied by 10 -3 . When the shading changes an extra factor 10 _1 should be inserted (e.g. 
for n q = 0.1357, fj(0) = 12.6(3) x 10~ 3 , / n (16) = 9.30(53) x 10" 4 , £(35) = 9.91(83) x 10" 5 , 
/i(68) = 9.0(1.4) x 10 -6 and so on). The results listed here correspond to the heavy-quark 
action, which we called in the text HYP-2 2 . Similar tables for other actions can be obtained 
from the authors. All results are given in lattice units. 
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